Low-CJ scaling, duality, and the EMC effect 
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High energy lepton scattering has been the primary tool for mapping out the quark distributions 
of nucleons and nuclei. Data on the proton and deuteron have shown that there is a fundamental 
connection between the low and high energy regimes, referred to as quark-hadron duality. We 
present the results of similar studies to more carefully examine scaling, duality, and in particular 
the EMC effect in nuclei. We extract nuclear modifications to the structure function in the resonance 
region, and for the first time demonstrate that nuclear effects in the resonance region are identical 
to those measured in deep inelastic scattering. With the improved precision of the data at large x, 
we for the first time observe that the large- a; crossover point appears to occur at lower x values in 
carbon than in iron or gold. 

PACS numbers: 25.30.Fj, 13.60.Hb 



I. INTRODUCTION 

Extensive measurements of inclusive lepton-nucleus 
scattering have been performed in deep inelastic scat- 
tering (DIS) kinematics. In DIS kinematics, where both 
the four-momentum transfer, Q, and the energy transfer, 
v, are sufficiently large, the extracted structure function 
exhibits scaling, i.e. is independent of Q 2 except for the 
well understood logarithmic QCD scaling violations. In 
this region, the structure function is interpreted as an in- 
coherent sum of quark distribution functions, describing 
the motion of the quarks within the target. 

Such measurements have unambiguously shown that 
the nuclear structure functions deviate from the pro- 
ton and neutron structure functions. Such modifications, 
termed the EMC effect after the first experiment to ob- 
serve them 0, demonstrate that the nuclear quark dis- 
tribution function is not just the sum of the proton and 
neutron quark distributions. Within two years of the 
first observation, hundreds of papers were published on 
the topic. After 20 years of experimental and theoretical 
investigation, the effect still remains a mystery. For de- 
tailed reviews of the data and models of the EMC effect, 
see Refs. 0, 01 

Existing measurements of the EMC effect indicate lit- 
tle Q 2 dependence, and an A dependence in the magni- 
tude, but not the overall form, of the structure function 
modification in nuclei. The nature of the modifications 
in nuclei depends primarily on Bjorken-x, x = Q 2 /2Mv, 
where in the parton model x is interpreted as the momen- 
tum fraction of the struck quark, and the nuclear effects 
are divided into four distinct regions. In the shadow- 
ing region, x < 0.1, the structure function is decreased 
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in nuclei relative to the expectation for free nucleons. 
In the anti-shadowing region, 0.1 < x < 0.3, the struc- 
ture function shows a small nuclear enhancement. For 
0.3 < x < 0.7, referred to as the EMC effect region, 
the nuclear structure function shows significant deple- 
tion. Finally, there is a dramatic enhancement as x in- 
creases further, resulting from the increased Fermi mo- 
tion of the nucleons in heavier nuclei. 

Explanations of the EMC effect are hampered by the 
lack of a single description that can account for the 
nuclear dependence of the quark distributions in all of 
these kinematic regimes. Here, we will limit ourselves 
to x > 0.3, the region where valence quarks dominate. 
Even in this limited region, there is not a single ex- 
planation that can completely account for the observed 
nuclear structure function modifications. If the nuclear 
structure function in this region is expressed as a convo- 
lution of proton and neutron structure functions, there 
are two alternative approaches used to describe the ob- 
served medium effect: (1) incorporating nuclear physics 
effects that modify the energy-momentum behavior of 
the bound proton with respect to the free proton, or 
(2) incorporating changes to the internal structure of the 
bound proton. It has been argued, most recently in Q, 
that the binding of nucleons alone can not explain the 
EMC effect. In addition, several attempts to explain the 
EMC effect in terms of explicit mesonic components ap- 
pear to be ruled out due to limits set by Drell-Yan mea- 
surements 0- Hence, the EMC effect may be best de- 
scribed in terms of modifications to the internal structure 
of the nucleon when in the nuclear environment. 

We note that while the EMC effect has been mapped 
out over a large range of x, Q 2 , and A, information is still 
rather limited in some regions. There are limited data 
on light nuclei (A < 9), and almost no data in the DIS 
regime at extremely large x, where the quark distribu- 
tions in nuclei are enhanced due to the effects of binding 
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and Fermi motion. Since binding and Fermi motion im- 
pact the EMC ratios for all x values, it is important to 
be able to constrain these effects in a region where other, 
more exotic, explanations are not expected to contribute. 
It should be possible to learn more about the EMC effect 
at large x by taking advantage of the extended scaling of 
structure functions in nuclei 0,0- In this paper, we at- 
tempt to quantify the deviations from perturbative scal- 
ing at large x, with the goal of improving measurements 
of the structure functions and the EMC ratios at large x. 



II. SCALING OF THE NUCLEAR STRUCTURE 
FUNCTION 

Inspired by a recent series of electron scattering exper- 
iments in Hall C at Jefferson Lab, we revisit the issues 
of scaling in nuclear structure functions and the EMC 
effect. The Hall C data are at lower invariant mass 
W, W 2 = Mp + 2M p v(l - x), and therefore higher x, 
than data thus far used to investigate the EMC effect. 
Most notably, these new data are in the resonance re- 
gion, W 2 < 4 GeV 2 . In the DIS region, W 2 > 4 GeV 2 , 
the Q 2 dependence of the structure functions is predicted 
by perturbative QCD (pQCD), while additional scaling 
violations, target mass corrections and higher twist ef- 
fects, occur at lower Q 2 and W 2 values. Thus, data in 
the resonance region would not naively be expected to 
manifest the same EMC effect as data in the deep inelas- 
tic scaling regime. The effect of the nuclear medium on 
resonance excitations seems non-trivial, and may involve 
much more than just the modification of quark distribu- 
tions observed in DIS scattering from nuclei. 

However, while resonance production may show differ- 
ent effects from the nuclear environment, there are also 
indications that there is a deeper connection between 
inclusive scattering in the resonance region and in the 
DIS limit. This connection has been a subject of inter- 
est for nearly three decades since quark-hadron duality 
ideas, which successfully described hadron-hadron scat- 
tering, were first extended to electroproduction. In the 
latter, Bloom and Gilman Q showed that it was possible 
to equate the proton resonance region structure function 
F%{v, Q 2 ) at low Q 2 to the DIS structure function F 2 {x) 
in the high-Q 2 scaling regime, where F% is simply the 
incoherent sum of the quark distribution functions. For 
electron-proton scattering, the resonance structure func- 
tions have been demonstrated to be equivalent on average 
to the DIS scaling strength for all of the spin averaged 
structure functions (F\, F?, F t,) 0, ^J, and for some 
spin dependent on es ( Ai ) (for a review of duality 
measurements, see |12j). 

The goal of this paper is to quantify quark-hadron du- 
ality in nuclear structure functions and to determine to 
what extent this can be utilized to access poorly under- 
stood kinematic regimes. While the measurements of du- 
ality from hydrogen indicate that the resonance structure 
function are on average equivalent to the DIS structure 
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FIG. 1: (Color online) The F2 structure function per nu- 
cleon vs £ for hydrogen (top), deuterium (middle), and 
iron(bottom). For the hydrogen and deuterium data (0.8 
< Q 2 < 3.3 GeV 2 ), the elastic (quasielastic) data have been 
removed. For the iron data (Q 2 < 5.0 GeV 2 ), a cut of 
W 2 > 1.2 GeV 2 is applied to remove the quasielastic peak. 
The curves are the MRST [3 (solid) and NMC (dashed) 
parameterizations of the structure functions at Q 2 = 4 GeV 2 , 
with a parameterization of the EMC effect 15] applied to 
produce the curve for iron. 



functions, it has been observed that in nuclei, this aver- 
aging is performed by the Fermi motion of the nucleons, 
and so the resonance region structure functions yield the 
DIS limit without any additional averaging |y, Q • 

FigureQ]shows the structure functions for hydrogen [j| , 
deuterium |^| , and iron Q , compared to structure func- 
tions from MRST |f| and NMC [3 parameterizations. 
Each set of symbols represents data in a different Q 2 
range, with the highest Q 2 curves covering the highest £ 
values. Note that the data are plotted as a function o f 
the Nachtmann variable, £ = 2x/(l + yjl + 4M 2 x 2 /Q 2 ), 
rather than x. In the limit of large Q 2 , £ x, and so 
£ can also be used to represent the quark momentum in 
the Bjorken limit. At finite Q 2 , the use of £ reduces scal- 
ing violations related to target mass corrections ^3] . The 
difference between £ and x is often ignored in high energy 
scattering or at low x, but cannot be ignored at large x or 
low Q 2 . The goal is to examine ^-scaling to look for any 
significant scaling violations beyond the known effects of 
perturbative evolution and target mass corrections. Ex- 
amining the scaling in terms of £ instead of x is only 
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FIG. 2: (Color online) F2 structure function per nucleon vs. 
Q 2 for deuterium at fixed values of £. Dashed lines show a 
logarithmic Q 2 dependence, with the value of d lni^/dlnQ 2 
determined at each £ value from SLAC data at high Q 2 (up to 
20 GeV 2 ). The solid lines denote W 2 =2.0 and 4.0 GeV 2 . The 
combined statistical and systematic uncertainties are shown. 
The hollow symbols are data from SLAC |Tg|. while the solid 
symbols are from Jefferson Lab 



an approximate way of applying target mass corrections, 
but it is a reasonable approximation to a more exact cor- 
rection p"?l | in the case of the proton, and the appropriate 
prescription for target mass corrections in nuclei is not 
as well defined. 

The transition from scaling on average in the proton 
to true scaling for nuclei is clearly visible. There is signif- 
icant resonance structure visible in hydrogen, but on av- 
erage the structure function reproduces the scaling curve 
to better than 2% globally and 5% locally around each 
resonance for Q 2 > 1 GeV 2 0- For deuterium, Fermi 
motion and other medium effects broaden the resonances 
to the point where only the A resonance has a clear peak, 
and the data at higher W 2 values, while still in the reso- 
nance region, is indistinguishable from the scaling curve 
except at the lowest Q 2 values. For the iron data, taken 
at somewhat higher Q 2 values, even the A is no longer 
prominent, and deviations from pQCD predictions are 
small, and limited to the tail of the quasielastic peak. 

We can study the quality of scaling in the resonance 
region more directly by examining the Q 2 dependence of 
the structure function at fixed £. Figure [21 shows the Q 2 
dependence of the structure function for several values of 
£. Above W 2 — 4 GeV 2 , the data are in the DIS region 
and the Q 2 dependence is consistent with the logarithmic 
Q 2 dependence from QCD evolution (dashed lines). Even 
at lower W 2 , where the data are in the resonance region, 
scaling violations are small. Above Q 2 = 3 GeV 2 , the 
data deviate from the logarithmic Q 2 dependence by <J 
10%, even down to W 2 = 2 GeV 2 . Data on heavier nuclei 
show a similar extended scaling in the resonance region, 
as seen in Fig. 3 of Ref . Q , although the largest scaling 
violations (in the vicinity of the quasielastic peak) are 



smaller, due to the increased Fermi smearing. 

Analyses of duality for the proton 0] and for nu- 
clei [13, El show that the moments of the structure func- 
tion, M n — J x n ~ 2 F2(x, Q 2 )dx is the nth moment, follow 
perturbative QCD evolution down to Q 2 ss 2 GeV 2 for 
the proton and to even lower values, Q 2 <, 1 GeV 2 , for 
nuclei. The fact that the moments follow the perturba- 
tive behavior is consistent with the observation that the 
structure function in Figs. ^ an d ID are i on average, in 
agreement with the perturbative structure function. 

The data indicate relatively small deviations from 
pQCD for Q 2 > 3 GeV 2 at all values of £ measured. 
These deviations decrease as Q 2 increases, making the 
nuclear structure functions at large £ consistent with the 
perturbative dependence even at values of W 2 well below 
the typically DIS limit. The limited kinematics coverage 
of the existing data make it difficult to precisely map out 
deviations from perturbative evolution. There is a large 
gap in Q 2 between the JLab data shown here and the 
SLAC measurements at large Q 2 . The situation will be 
improved by the recently completed measurements from 
JLab experiments E03-103 and E00-116 [HIH, which 
will provide more complete £ coverage over a wide range 
in Q 2 . In the meantime, Fig.|5]indicates that for Q 2 <; 3- 
4 GeV 2 , one can relax the usual DIS requirement that 
W 2 > 4 GeV 2 , and the structure functions measured in 
the resonance region will still provide a good approxima- 
tion to the DIS structure functions. 

Even with the uncertainties arising from possible 
higher twist contributions, data at large £ can signifi- 
cantly improve our knowledge of the high-£ nuclear struc- 
ture functions. There is very little DIS data for £ ;> 0.8, 
and no existing facility has the combination of energy 
and luminosity necessary to make precise measurements 
of in this regime, which requires Q 2 > 15(36) GeV 2 for 
£ > 0.8(0.9). If we can set reasonable limits on scaling 
violations at fixed £ due to possible higher twist contri- 
butions, we can provide useful data in this region where 
the few existing measurements of the EMC effect have 
10-20% uncertainties. 



III. STRUCTURE FUNCTION RATIOS 

Because of the difficulty in making precise measure- 
ments in the DIS region at large £, existing measurements 
of the EMC effect in this region are very poor. At large 
£, the EMC effect should be dominated by binding ef- 
fects and Fermi motion. Constraining these effects will 
allow a better separation of these "conventional" nuclear 
effects, which are important at all £ values, from more 
exotic effects that have been used to explain the EMC 
effect at lower £. 

We can examine the EMC effect in the resonance re- 
gion using recent measurements of inclusive scattering 
from deuterium, carbon, iron, and gold. For these data, 
we take the cross section ratio of iron to deuterium in the 
resonance region for the highest Q 2 measured (Q 2 ~ 4 
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GeV 2 ), requiring W 2 > 1.2 GeV 2 to exclude the region 
very close to the quasielastic peak. 

There are small differences between the analyses of 
the SLAC and JLab data which had to be addressed 
to make a precise comparison. First, the SLAC and 
BCDMS ratios were extracted as a function of x rather 
than £. Because the conversion from x to £ depends on 
Q 2 , we can only compare ratios extracted at fixed Q 2 
values. Thus, for E139 we use the "coarse-binned" ra- 
tios, evaluated at fixed Q 2 , rather than "fine" x binning, 
which were averaged over the full Q 2 range of the exper- 
iment. Coulomb corrections were applied in the analy- 
sis of the JLab data [24]], but not the SLAC data. The 
SLAC data shown here include Coulomb corrections, de- 
termined by applying an offset to the incoming and out- 
going electron energy at the reaction vertex [24}, due to 
the Coulomb field of the nucleus. The correction fac- 
tor is <0.5% for carbon, and (1.5-2.5)% for gold. The 
JLab and SLAC ratios are corrected for neutron excess, 
assuming a n ja p = (1 — 0.8£). 

Figure [3] shows the cross section ratio of heavy nuclei 
to deuterium for the previous SLAC E139 [H, E87 
and BCDMS |2q DIS measurements, and for the JLab 
E89-008 0, l24| data in the resonance region. The size 
and £ dependence of nuclear modifications in the JLab 
data agrees with the higher Q 2 , W 2 data for all targets. 
Table |U shows the ratios extracted from the JLab data. 

The agreement of the resonance region data with the 
DIS measurement of the EMC effect, which directly mea- 
sures the modification of quark distributions in nuclei, is 
quite striking. There is no a priori reason to expect that 
the nuclear effects in resonance production would be sim- 
ilar to the effects in scattering from quarks. However, it 
can be viewed as a natural consequence of the quantita- 
tive success of quark-hadron duality 0, • As seen m 
Fig. 2] the structure functions for nuclei show little devi- 
ation from pQCD, except in the region of the quasielastic 
peak (and A resonance at low Q 2 ). As Q 2 increases, the 
deviations from pQCD decrease as quasielastic scattering 
contributes a smaller fraction of the cross section. In ret- 
rospect, given the lack of significant higher twist contri- 
butions, combined with the fact that any A-independcnt 
scaling violations will cancel in the ratio, it is perhaps not 
surprising that the resonance EMC ratios are in agree- 
ment with the DIS measurements. 

While it is difficult to precisely quantify the higher 
twist contributions with the present data, we can esti- 
mate their effect by looking at low W 2 and Q 2 , where 
the higher twist contributions are much larger. At Q 2 f=a 
2 GeV 2 and W 2 « M\ : the scaling violations (beyond 
target mass corrections) for deuterium are as large as 
50%, as seen in Fig. ^ However, if one takes the iron 
and deuterium data from Ref. , averages the structure 
function over the A region and then forms the EMC ra- 
tio, the result differs from the ratio in the DIS region by 
less than 10%. The decrease in the effect of higher twist 
contributions is a combination of the fact that the con- 
tribution are reduced when averaged over an adequate 
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FIG. 3: (Color online) Ratio of nuclear to deuterium cross 
section per nucleon, corrected for neutron excess. The solid 
circles are Jefferson lab data taken in the resonance region 
(1.2 < W 2 < 3.0 GeV 2 , Q 2 w 4 GeV 2 ). The hollow diamonds 
are SLAC E139 data, the crosses are the SLAC E87 data, and 
the hollow squares are BCDMS data, all in the DIS region. 
The scale uncertainties for the SLAC (left) and JLab (right) 
data are shown in the figure. The curves show an updated 
version |27j | of the calculations from Ref. |2S| . 



region in W 2 |3,[T3|; and cancellation between the higher 
twist contributions in deuterium and iron. The same 
procedure yields 2-3% deviations from the EMC ratio 
if one looks in the region of the Sn or Pi 5 resonances, 
where the scaling violations in the individual structure 
functions are smaller to begin with. 

For the ratios in Fig. El we expect even smaller higher 
twist effects because the data is nearly a factor of two 
higher in Q 2 and is above the A except for the very 
highest £ points. At higher Q 2 , the higher twist con- 
tributions in the individual structure functions become 
smaller, while averaging over the resonance region be- 
comes less important as the resonances become less 
prominent. Thus, we expect that higher twist contri- 
butions for these data will be smaller than the the 2-3% 
effect (<10% near the A) observed on the EMC ratio at 
Q 2 « 2 GeV 2 . If so, the higher twist corrections will 
be small or negligible compared to the large statistical 
uncertainty in previous measurements, and this data can 
be used to improve our knowledge of the EMC effect at 
large £. 
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W 2 


(ac /(TD)ia 


(f"Fe/o"D)is 


{(JAilI <JD)is 




GeV 2 








0.592 


2.86 


0.921±0.012 


0.901±0.013 


0.844±0.013 


0.613 


2.70 


0.926±0.013 


0.905±0.014 


0.871±0.015 


0.633 


2.55 


0.882±0.013 


0.872±0.015 


0.843±0.016 


0.654 


2.39 


0.908±0.014 


0.921±0.017 


0.865±0.017 


0.676 


2.22 


0.920±0.012 


0.881±0.012 


0.872±0.014 


0.697 


2.07 


0.948±0.011 


0.883±0.010 


0.874±0.013 


0.719 


1.91 


0.971±0.013 


0.917±0.012 


0.903±0.016 


0.741 


1.75 


1.000±0.016 


0.976±0.015 


0.938±0.020 


0.763 


1.59 


1.031±0.019 


0.964±0.017 


0.937±0.023 


0.786 


1.43 


1.139±0.022 


1.013±0.020 


1.043±0.028 


0.810 


1.26 


1.195±0.014 


1.133±0.015 


1.138±0.024 



TABLE I: Isoscalar EMC ratios for carbon, iron, and gold in 
the resonance region extracted from the data of Ref. • W 2 
is calculated using the nucleon mass, rather than the nuclear 
mass. 



IV. THE EMC EFFECT AT LARGE x (£) 

A careful examination of the crossover point at large £, 
where the ratio (oa/cd)^ becomes larger than unity, re- 
veals that this appears to occurs at larger £ for heavy nu- 
clei than for light nuclei. This behavior is consistent with 
the SLAC data, but the large-£ coverage of the previous 
measurements was insufficient to make a clear statement 
about the crossover point. This observation contradicts 
the argument that the dramatic enhancement at large £ 
is simply due to increased Fermi motion in heavy nuclei 
relative to deuterium. In this simple picture, the slight 
increase in Fermi motion as one goes to heavier nuclei 
would lead to an earlier onset of this enhancement. While 
the Coulomb corrections and neutron excess corrections 
do affect the A dependence, the uncertainties on these 
corrections are not large enough to explain the observed 
differences between carbon and heavier nuclei. 

Within the convolution formula of proton and neutron 
structure functions, this crossover comes about due to 
counteracting contributions at large £ of the av erag e nu- 
cleon binding energy and average kinetic energy |29| , and 
is hardly expected to change for A > 10. More detailed 
calculations have been done to determine the effect of 
binding and Fermi motion, but except for calculations of 
3 ' 4 He, most of these calculations were performed for a 
Fermi gas model or for infinite nuclear matter, with the 
density varied to approximate the finite nuclei. These 
calculations do not show any significant change in the £ 
dependence for heavy (A > 10) nuclei. Because of the 
lack of precise data at large x, especially with respect to 
the ^-dependence, realistic models of the nuclear struc- 
ture were generally not considered to be necessary. 

The effect we observe was predicted in a calculation by 
Gross and Liuti using a manifestly covariant form of 
the convolution formula. The most significant difference 



was an additional binding correction due to the explicit 
dependence of the bound nucleon structure function on 
the momentum of the bound nucleon. Their calculation 
predicts a shift in the high-£ crossover point between car- 
bon and iron, somewhat larger than is observed in the 
data. Another calculation including A-dependent nu- 
clear spectral functions [3(j also gave an A-dependent 
crossover point at large £. However, this calculation had 
the crossover point moving to lower £ values for heavier 
nuclei. 

Other models of the EMC effect have looked at physics 
beyond Fermi motion and binding. As with the bind- 
ing models, they generally did not attempt to reproduce 
the detailed A-dependence, and instead evaluated the 
EMC effect for nuclear matter as a function of density. 
Most of these models were designed to describe the excess 
strength at lower £, and in general they do not signifi- 
cantly impact the structure function at large £. Thus, the 
addition of improved EMC ratio measurements at large 
£ and the observation of an A dependence to the high-£ 
behavior is most important in constraining the portion 
of the EMC effect that is related to binding. One can 
see from the calculations shown in Fig. that the effects 
of binding and Fermi motion are important over the en- 
tire £ region, and not just at the largest £ values. These 
conventional nuclear effects must be well constrained to 
establish a reliable baseline before one can isolate any 
additional nuclear modification at lower x values that 
might require a more exotic explanation. Improved data 
at large £ and for a variety of nuclei should allow for tests 
of the prescriptions chosen for binding and Fermi motion, 
and thus provide a more reliable baseline for models of 
the EMC effect. 

The modified x dependence in carbon also appears to 
contradict the conclusions of a recent effective field the- 
ory calculation of the EMC effect that predicts factoriza- 
tion of the A and x dependence, and thus the universal- 
ity of the x dependence [3j]|. The change in the high- a: 
crossover in the present data is small, but a recent mea- 
surement of the EMC effect for 3 He and 4 He ^ will 
provide a more sensitive test of the universality of the x 
dependence. 



V. CONCLUSIONS 

This analysis provides not only an increase in the £ 
range of the measurements of nuclear structure functions, 
but also the first observation of an A dependence of the 
high-£ behavior of the EMC effect. Measurements uti- 
lizing higher energy beams will extend measurements of 
the EMC effect to even larger £ values. Based on the 
results shown here, the uncertainties on extracting the 
EMC effect at large £ due to higher twist contributions 
will be small, if not negligible, compared to the uncertain- 
ties of existing data. A recent measurement at Jefferson 
Lab [2^1 will extend measurements of the EMC effect 
to larger £ values and to light nuclei, where few-body 
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calculations can be performed with significantly smaller 
uncertainties coming from uncertainties in the nuclear 
structure. The calculations of Ref. predict that the 
high-i; crossover point in carbon occurs at lower £ than 
in heavier nuclei, as was observed in the data. However, 
they predict a crossover at much larger £ for He. Sim- 
ilarly, Ref. also predicts a different high-£ behavior 
in He than in heavy nuclei, and in addition predicts a 
significant difference between 3 He and 4 He. 

Similar investigations of duality and scaling in po- 
larized and separated structure functions are under- 
way [33L I34L |33 • If duality in these processes is quan- 
titatively as successful as in this case, this will have a 
similar impact on our ability to measure high-£ polarized 
structure functions. 

In conclusion, we present the first extraction of the nu- 
clear dependence of the inclusive structure function in the 
resonance region. The data are in agreement with previ- 
ous measurements of the nuclear dependence of the quark 
distributions in DIS scattering measurements of the EMC 
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